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ABSTRACT 

Context. The Great Nebula in Carina is a giant H II region and a superb location in which to study the physics of violent massive star 
formation, but the population of the young low-mass stars remained very poorly studied until recently. 

Aims. Our aim was to produce a near-infrared survey that is deep enough to detect the full low-mass stellar population (i.e. down to 
a; 0. 1 M and for extinctions up to Ay as 15 mag) and wide enough to cover all important parts of the Carina Nebula complex (CNC), 
including the clusters Tr 14, 15, and 16 as well as the South Pillars region. 

Methods. We used HAWK-I at the ESO VLT to survey of the central as 0.36 deg 2 area of the Carina Nebula. These data reveal more 
than 600000 individual infrared sources down to magnitudes as faint as J « 23, H « 22, and K s « 21. The results of a recent deep 
X-ray survey (which is complete down to stellar masses of ~ 0.5 - 1 M G ) are used to distinguish between young stars in Carina and 
background contaminants. We analyze color-magnitude diagrams (CMDs) to derive information about the ages and masses of the 
low-mass stars. 

Results. The ages of the low-mass stars agree with previous age estimates for the massive stars. The CMD suggests that « 3200 
of the X-ray selected stars have masses of M, > 1 M ; this number is in good agreement with extrapolations of the field IMF based 
on the number of high-mass (M, > 20 M G ) stars and shows that there is no deficit of low-mass stars in the CNC. The HAWK-I 
images confirm that about 50% of all young stars in Carina are in a widely distributed, non-clustered spatial configuration. Narrow- 
band images reveal six molecular hydrogen emission objects (MHOs) that trace jets from embedded protostars. However, none of the 
optical HH objects shows molecular hydrogen emission, suggesting that the jet-driving protostars are located very close to the edges 
of the globules in which they are embedded. 

Conclusions. The near-infrared excess fractions for the stellar population in Carina are lower than typical for other, less massive 
clusters of similar age, suggesting that the process of circumstellar disk dispersal proceeds on a faster timescale in the CNC than in 
the more quiescent regions, most likely due to the very high level of massive star feedback in the CNC. The location of all but one of 
the known jet-driving protostars at the edges of the globules adds strong support to the scenario that their formation was triggered by 
the advancing ionization fronts. 

Key words. Stars: formation - Stars: mass function - Stars: circumstellar matter - Stars: pre-main sequence - ISM: individual 
objects: NGC 3372 - open clusters and associations: individual: Tr 14, Tr 15, Tr 16 

1 . Introduction type stars profoundly influence their environment by their strong 

ionizing radiation, powerful stellar winds, and, finally, by su- 



H ■ Most stars in the Galaxy are thought to be born in massive pernova eX pi osio ns. This fee dback can disperse the surrounding 

. star-forming regions (e.g., | Blaauw| | 1964t | Miller & Scalo||1978t natal molecular clouds (e.g. JFrever et allMl) . and thus termi- 

i Bnceno et al. | |2007|), and therefore in close proximity to mas- nate the star formation process ( = negative feedback). However, 

sive stars. This also applies to the origin of our solar system, ion i Z ation fronts and expanding superbubbles can also compress 

for which recent investigations found convincing evidence that nearby clouds and thereby trigger the f orm at ion of new gen- 

it formed in a large cluster, consisting of (at least) several thou- eratio ns of stars (= positive feedback; e.g., iGritschneder etall 

sand stars, and that the original solar nebula was directly affected ffififr | Pre ibisch & Zinneckeril2007r l . These processes determine 

by nearby massive stars (e.g., | Adams || 2010| ). As a consequence, me key outputs from star formation , including the stellar mass 

the role of environment is now an essential topic in studies of fun ction, the total star formation efficiency, as well as th e evolu- 

star and planet formation. The presence of high-mass stars can tion of circumste ii ar disks around the young stars (e.g., Iciarkel 

lead to physical conditions that are vastly different from those in ffiM and the frequency of planetary formation. While this gen- 

regions where only low-mass stars form. The very luminous O- eral picture is now well established, the details of the feedback 

processes, i.e. cloud dispersal on the one hand, and triggering 

* Based on observations collected with the HAWK-I instrument at the of star formation on the other hand, are still only poorly under- 
VLT at Paranal Observatory, Chile, under ESO program 60.A-9284(K). 
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stood. A fundamental problem is that nearly all massive star- 
forming clusters with high levels of feedback are quite far away 
and therefore difficult to study. At distances greater than 5 kpc, 
the detection and characterization of the full stellar populations 
is very difficult (if not impossible); often, only the bright highl- 
and intermediate-mass stars can be studied, leaving the low-mass 
(M < 1 Mq) stars unexplored or even undetected. However, 
since the low-mass stars constitute the vast majority of the stellar 
population, a good knowledge of the low-mass stellar content is 
essential for any inferences on the initial mass function and for 
understanding the nature of the star formation process. Detailed 
information about the low-mass stars and their protoplanetary 
disks is also crucial for investigations of the interaction between 
the high- and low-mass stars. The intense UV radiation from 
massive stars may remove considerable amounts of the circum- 
stellar material from nearby young stellar objects (YSOs, here- 
after). This may limit the final masses of the low-mass stars (see 
IWhitworth & Zinne cker 200 % and should also affect the forma- 
tion of planets (see, e.g., Throop & Ballvll2"b05h . 

In thi s context, the Great N ebula in Carina (NGC 3372; 
see, e.g., ISmifh & Brook s 2008, for an overview) provides a 
unique target for studies of massive star feedback. At a very 
well known and moderate distance of 2.3 kpc, the Carina Nebula 
Complex (CNC, hereafter) represents the nearest southern re- 
gio n with a larg e massive stellar population (65 O-type stars; 
see ISmithl 120061) . Among these are several of the most mas- 
sive (M ^ 100 M Q ) and luminous stars known in our Galaxy, 
e.g., the famous Luminous Blue Variable 77 Car, the 02 If* star 
HD 93129Aa, several 03 main sequence stars and Wolf-Rayet 
stars. The CNC is the most nearby region that samples the 
top of the stellar mass function. The presence of stars with 
M k, 100M Q implies that the level of feedback in the CNC is 
already close to that in more extreme extragalactic starburst re- 
gions, while at the same time its comparatively moderate dis- 
tance guarantees that we still can study details of the cluster and 
cloud structure at good enough spatial resolution and detect and 
characterize the low-mass stellar populations. Due to this unique 
combination of properties, the CNC represents the best galactic 
analog of giant extragalactic H II and starburst regions. 

Most of the very massive stars in the CNC reside in sev- 
eral loose clusters, including Tr 14, 15, and 16, for which have 
ages between ~ 2 an d ~ 8 Myr have been estimated (see 
ISmifh & Brooksl 120081) . The CNC is thus often denoted as a 
"cluster of clusters". In the central region, the molecular clouds 
have already been largely dispersed by the feedback from the 
numerous massive stars. Southeast of 77 Car, in the so-called 
"South Pillars" region, the clouds are eroded and shaped by the 
radiation and winds from 77 Car and Tr 16, giving rise to nu- 
merous giant dust pillars, which feature very prominently in the 
mid-infrared image s made with the Spitzer Space Observatory 
dSmith et al.l2010b|). The complex contains more than 10 5 M Q o f 
gas and dust (see Yonekura et al.ll2005l iPreibisch et al.ll201 1 ah . 
and deep infrared observations show clear evidence of ongo- 
ing star fo rmation in these clo uds. Several deeply embedded 
YSOs (e.g.. lMottram et al.l2007l) a nd a spectacular yo ung cluster 
(the "Treasure Chest Cluster"; see ISmith et aDl2005l) have been 
found. A deep HST Ha imaging s urvey revealed dozens of jet- 
driving YSOs (Smith et a"Dl2010al) . and Spitzer surveys located 
numerous embedded intermediate-mass protostars throughou t 
the Carina complex (ISmithet al.l 1201 Obt iPovich et all 1201 lab . 
The formation of this substantial population of very young (< 
1 Myr), partly embedded stars was probably triggered by the ad- 
vancing ionization fronts that originate from the (several Myr 
old) high-mass stars. The CNC provides an excellent target to 



1998). 



investigate these effects on the formation and evolution of low- 
mass stars (and their forming planetary systems) during the first 
few Myr. 

While the un-obscured population of high-mass stars (M > 
20 Mq) in the CNC is well known and characterized, the (much 
fainter) low-mass (M < 1 M Q ) stellar population remained 
largely unexplored until now. One reason for this is the faint- 
ness of the low-mass stars: at a distance of 2.3 kpc, an age 
of 3 Myr, and assumin g typical extinctions of Ay = 3.5 mag 
(IPreibisch et al.ll201 lbl) . stars with masses of 0.5 M [0.1 M Q ] 
are predicted to have magnitudes of V - 23.6 [26.6 
17.0 [19.2], and K = 15.5 [17.7] (see iBaraffe etal] 
The existing near-infrared (NIR) observations were either 
too shallow to detec t the full low-mass population (e.g., 
ISanchawala et al.ll2007l). or covere d only very small parts of the 
complex (e.g. Ascenso et a l. 2007). 

This was the motivation for the deep wide-field NIR sur- 
vey presented in this paper. Our survey is > 5 mag deeper than 
the 2MASS data, deep enough to detect essentially all young 
stars in the CNC, and can be used to derive information about 
the ages and masses of the young low-mass stars. A funda- 
mental restriction is that the NIR data alone do not allow us to 
distinguish young stars in the CNC from the numerous galac- 
tic field stars in the area (see discussion in Sect. 2). For this 
purpose, we use the recent results from the Chandra Carina 
Complex Project (CCCP), that has mapped the entire extent of 
the CNC (~ 1.4deg 2 ) in X-rays. A comple te overview of the 
CCCP can be found in lTownslev et alJ d201 lh and associated pa- 
pers in a Special Issue of the Astrophysical Journal Supplements. 
The CCCP data led to the detection of 14 368 individual X-ray 
sources. Based on a statistical analysi s of the properties of each 
individual source. iBroos et al" classified 10714 of these 

as young stars in the CNC. The identification of HAWK-I in- 
frared counterparts of the X-ray sources in the CNC and a few 
basic aspects of the X-ray sele cted population of you ng stars 
have been recently presented in IPreibisch et al.l (1201 lbl) . In this 
paper we discuss the complete set of HAWK-I data (Sec. 2) and 
present a detailed analysis of the NIR color-magnitude diagrams 
(Sec. 3). In Section 4 we analyze the NIR excess fractions for the 
different clusters in the CNC. In Section 5 we draw global impli- 
cations on the size of the stellar population and the IMF of the 
CNC. Section 6 contains the results of our search for protostellar 
jets in H2 narrow-band images. In Sec. 7 we look at the spatial 
distribution of the young stars and list newly discovered stellar 
clusters in the CNC. Global conclusions on the star formation 
process in the CNC are discussed in Sec. [8] 

2. HAWK-I observations and data analysis 

The NIR imager HAWK-I (see lKissler-Patig et alJl20"08T) at the 
ESO 8 m Very Large Telescope is equipped with a mosaic of four 
Hawaii 2RG 2048 x 2048 pixel detectors with a scale of 0.106" 
per pixel. The camera has a field of view on the sky of 7.5' x 
7.5' with a small cross-shaped gap of ~ 15" between the four 
detectors. The observations of the CNC were performed from 
24 to 3 1 January 2008 in service mode as part of the scientific 
verification program for HAWK-I. In addition to the standard 
broadband J, H, and K s filters, we used the H 2 (2. 109-2. 139 pm) 
and Bry (2.150 - 2.181 pm) narrow-band filters. 

Our survey consists of a mosaic of 24 contiguous HAWK-I 
fields, covering a total area of about 1280 square-arcminutes; it 
includes the central part of the Nebula with 77 Car and Tr 16, the 
clusters Tr 14 and Tr 15, and large parts of the South Pillars (see 

Fig.iD. 
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Fig. 1. Negative gray scale representation of an optical image of the Carina Nebula taken from the red Digitized Sky Survey plates, 
with superposed outline of the HAWK-I mosaic pattern. The individual HAWK-I mosaic fields are shown as boxes with a size of 
7.5' X 7.5'. The clusters Tr 14, Tr 15, Tr 16, Cr 232, Bo 11, and the Treasure Chest Cluster (TCC) are marked. The 2' diameter blue 
circles mark the positions of the four embedded clusterings discussed in Sect. [6] and the six small red circles indicate the positions 
of the six Molecular Hydrogen Emission Objects discussed in 17. II North is up and east to the left. A grid of J2000 coordinates is 
also shown. 



The total observing times for each mosaic position were 12, 
8, and 5 minutes in J, H, and K s . Additional 5 minutes were 
spent with each of the H2 and Bry narrow band filters. In order 
to account for the gaps of the detector array, we used a 5-point 
dither pattern with offsets of (0,0), (-40", +40"), (-40", -40"), 
(+40", -40"), and (+40", +40") at each individual mosaic po- 
sition. Details of the observing parameters are listed in Table Q] 

Thanks to the very good seeing conditions during the ob- 
servations, most of the HAWK-I images are of a very high im- 
age quality. The FWHM of the Point-Spread-Function in the 
HAWK-I images is typically 0.6" - 0.8". This good image qual- 
ity led to the detection of interestin g structures such as the edge- 
on circumstellar disk described in ( Preibisch et al.ll201 lch . 



Table 1. Observing parameters 



Filter 


Offsets 


Integrations 


Integration 
times [s] 


Total exposure 
time [s] 


J 


2x5 


24 


3 


720 


H 


2x5 


16 


3 


480 


K s 


5 


20 


3 


300 


H 2 


5 


2 


30 


300 


Bry 


5 


2 


30 


300 



As examples for the quality of the HAWK-I data, we show 
comparisons of the HAWK-I and 2MASS images of selected re- 
gions in the CNC in Fig. [2] 
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Fig. 2. Comparison of HAWK-I (left) and 2MASS (right) //-band images of Bochum 11 (above) and the Treasure Chest cluster 
(below). 



2. 1 . Photometry and construction of the HAWK-I source 
catalog 

The procedures performed to derive photometry from the 
HAWK-I images and the construction of the source catalog is 
described in detail in iPreibisch et al.l (1201 lbl) . Here we briefly 
summarize the main points. All HAWK-I data were processed 
and calibrated by the Cambridge As tronomical Survey Unit us- 
ing pipeline software as described in llrwin et al. Id2004b. A pho- 
tometric analysis was performed for all images obtained in the 
three broad-band J, H, and K s , follow i ng clo sely the procedures 
described in detail in iHodgkin et al.l d2009l) . The photometric 
calibrators are drawn from stars in the 2MASS Point Source 
Catalog, which are present in large numbers (several hundred) 
in each HAWK-I mosaic field. The source lists of the individ- 
ual mosaic positions were then combined into total catalogs (per 
band) by removing duplicates from the overlap regions; this was 
achieved by keeping the source with the highest signal-to-noise 
ratio. The final catalog was generated by merging the source cat- 
alogs for the individual bands into a combined catalog. The max- 
imum allowed spatial offset for inter-band identifications of in- 
dividual objects was set to < 0.3". 

The final HAWK-I photometric catalog contains 600336 in- 
dividual objects. Most (502714) catalog objects are simultane- 
ously detected in the J-, H-, and the A^-band. Objects as faint as 
J ~ 23, H ~ 22, and K s ~ 21 are detected with S/N > 3. Typical 



values for the completeness limit across the field are 7 comp i ~ 21, 
#compi ~ 20, and A'.t.compi ~ 19; nearly all objects above these 
limits are S/N > 10 detections. Our survey is thus about 5 mag 
deeper than 2MASS and represents the largest and deepest NIR 
survey of the CNC obtained so far. 

The absolute photometric accuracy of the catalog was char- 
acterized by a comparison of the HAWK-I and 2MASS photom- 
etry for a sample of more than 700 stars. The standard devia- 
tions between HAWK-I and 2MASS magnitudes and colors are 
found to be cry = 0.11 mag, cr H = 0.08 mag, o~k s =0.11 mag, 
Cj-h = 0.12 mag, and cth-k s = 0.11 mag. The accuracy is 
mainly limited by the photometric quality of the 2MASS stars 
suitable for comparison and the strong and highly variable dif- 
fuse nebulosity in the area of the HAWK-I mosaic. 

2.2. Background contamination in the HAWK-I catalog 

Due to the CNC's location almost exactly on the galactic plane, 
the degree of field star contamination is very high. This causes a 
fundamental limitation of the usefulness of the infrared source 
catalog. The degree of contaminatio n was recen tl y quan tified 
in two Spitz e r studi es of the CNC. ISmifh et all (1201 Obi) and 
iPovich et all (l2011al) found that only * 3% of the ~ 50000 
Spitzer infrared sources in the CCCP survey area can be clas- 
sified as candidate YSOs. The contamination rate of our deep 
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Fig. 3. Color-Magnitude Diagram of all HAWK-I sources with 
NIR excesses. The solid lines show isochrones for ages of 1 Myr 
(red), 3 Myr (gree n), 10 Myr (blue), an d the ZAMS composed 
from the models of Baraffe et al for the mass range 0.02 

to 0.5 M and ISiess et al.l (120001) for the mass range 0.5 to 
7 M Q . The large solid dots mark the positions of 1 M Q stars, 
the asterisks those of 0.075 M e objects on the 1, 3, and 10 Myr 
isochrones. extinction vector for Ay = 10 mag. The arrows in- 
dicate reddening vectors for Ay = 10 mag and were com puted 
for an extinction law with Ry — 4 (see lPovic h et al. 201 lb). The 
dashed green line shows the 3 Myr isochrone reddened by that 
amount. Note that all bright objects with J < 12.5 are in the 
saturation/non-linearity regime of the HAWK-I data and their 
photometry is therefore unreliable. 



NIR HAWK-I data is presumably similarly high. This suggests 
that the vast majority (> 90%) of the ~ 600000 detected in- 
frared sources in our HAWK-I catalog must be foreground or 
background sources, unrelated to the CNC. 

2.3. Contamination and incompleteness of NIR excess 
selected samples 

This high degree of contamination requires a very efficient and 
reliable selection method to identify YSOs among the very nu- 
merous contaminants. A widely used method is to select objects 
with infrared excess emission, which is a tracer of circumstel- 
lar material, as YSO candidates. In order to test whether such 
an approach could be useful for our data, we analyzed the J-H 
vs. H—K s color-color diagram of the HAWK-I sources. Objects 
are classified as NIR excess sources if they lie at least > Icr^ot 
and more than 0.05 mag to the right and belo w the reddening 
band, based on the intrinsic colors of dwarfs (IBessel & Brettl 
[l988h and an extinction vector with slope 1 .73, and above J-H = 
0. In this way, 50950 of the 502714 HAWK-I sources with 
complete J, H, and K s photometry are classified as NIR excess 
sources. The color-magnitude diagram (CMD) of these NIR ex- 
cess sources is shown in Fig. [3] 



However, there are two reasons why we think that this 
excess-selected sample is not a useful sample of the YSOs in the 
CNC. First, due to the very high level of background contami- 
nation, even an excess selected sample will still be highly pol- 
luted by background objects. Several recent spectroscopic stud- 
ies of infrared sources in other star forming regions have shown 
that excess-selected YSO candidate samples can be strongly 
contaminated by background sources such as evolved Be stars, 
carbon stars, or planetary nebulae, which often show NIR col- 
ors very similar to those of YSOs (e .g., iMentuch et al.l 120091 
lOliveira et aLll2009t iRebull et al.llIoToh . Our CMD suggests this 
to be the case here too: a large majority (78%) of the NIR excess 
objects is found at very faint magnitudes (J > 19). Comparing 
the location of the sources in the CMD to pre-main sequence 
models shows that 61% of all NIR excess objects are below the 
expected reddening vector for 3 Myr old 0.075 M Q objects in 
Carina. While some of these faint objects may be deeply em- 
bedded YSOs, stars with edge-on circumstellar disks, or disk- 
bearing young brown dwarfs, the majority are very likely not 
young stars in the CNC but background objects. Furthermore, we 
note that the HAWK-I images are deep enough to detect numer- 
ous extragalactic objects. Since star-forming galaxies and AGN 
often display NIR colors similar to those of YSOs, this leads to 
even higher contamination rates in NIR excess selected samples 
of YSO candidate^ 

The second problem of a NIR excess selected YSO candidate 
sample is its incompleteness. It is well known that NIR excess 
emission in young stars disap pears on timescales of just a few 
Myr (e.g.. lBriceno et al.120071) : at an age of ~ 3 Myr, only ~ 50% 
of the young stars still show NIR excesses, and by ~ 5 Myr this is 
reduced to ~ 15%. Since the expected ages of most young stars 
in the CNC are several Myr, any excess-selected YSO sample 
will be highly incomplete. 

Considering these problems, we conclude that due to the 
very high degree of background contamination and the high ex- 
pected level of incompleteness of NIR excess selected YSO can- 
didate samples, the NIR data alone cannot provide a reasonably 
clean and complete sample of YSOs. Without additional infor- 
mation, the HAWK-I data alone can yield only very limited in- 
sight into the young stellar populations in the CNC. 



2.4. Combination of infrared and X-ray data 

Sensitive X-ray observations can provide a very good solution 
for the problem of identifying the young stars in an extended 
complex such as the CNC. X-ray sur veys detect the young stars 
by their strong X-ray emission (e.g., iFeigelson et al.ll2007l) and 
efficiently discriminate them from the numerous older field 
stars in the survey area. X-rays are equally sensitive to young 
stars which have already dispersed their circumstellar disks, 
thus avoiding the bias introduced when selecting samples based 
only on infrared excess. Many X-ray studies of star forming 
regions have demonstrated the s uccess o f this method (see, e.g. 
Preibisch. Zinnecker & Herbid [T9961 iPreibisch & Zinneckei 



2002; iBroos et al.ll2007l: iForbrich & Preibischll2007l IWang et al. 
l2010h . Also, the relations between the X-ray properties and 
basic stellar properties in young stellar populations are now 



1 Based on the number counts of star-forming g alaxies in the data 
from the UKIRT Infrared Deep Sky Survey data by lLane et all d2007f) 
we can expect the presence of about 4700 star-forming galaxies with 
K s < 20 in the HAWK-I field. From a purely statistical point-of-view, 
all of the 2157 NIR excess objects with K s > 20 in the HAWK-I catalog 
may well be extragalactic sources. 
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very well established from very deep X-ray observations such as 
the Chandra Orion U ltradeep Project (see iGetman et al.ll2.Q05t 
iPreibisch et"aT1l2005l) . 

The combination of the HAWK-I infrared source cata- 
log with the list of X-ray sources detected was an essen- 
tial asj)e£t^_of_the C/zant/ra Carina Complex Project (CCCP; 
see iTownslev et all 1201 lb . The d etails of the inf r ared - X-ray 
source matching is described in IPreibisch et al.l (1201 lbl) . and 
here we just briefly summarize the main aspects: The area of the 
HAWK-I mosaic covers 27% of the CCCP survey area and con- 
tains 52% of all 14 368 detected X-ray sources. The HAWK-I in- 
frared catalog provided infrared matches for 6583 of these X-ray 
sources. After adding matches from the 2MASS catalog for very 
bright stars that are saturated in the HAWK-I images, infrared 
counterparts were established for 6636 X-ray sources. 6241 of 
these have valid photometry in all three of the J-, H-, and K s - 
bands, and 6173 (9 3.0%) are classified as Carina members (see 
iBroos et~aT ] l2011bh . 

Using the NIR photometry of the X-ray selected Carina 
members, it was found that the typical extinction values for the 
diskless stars range from Ay -1.6 mag to ~ 6.2 mag (central 
80% percentile). While this clearly shows a considerable range 
of differential extinction between individual stars in the com- 
plex, it also implies that the vast majority of the X-ray detected 
young stellar population in the CNC shows rather moderate ex- 
tinction, Ay < 10 mag. The X-ray selected Carina members can 
thus be considered as a "lightly obscured" population of young 
stars. As will be discussed in Sec. [8] deep Spitzer data suggest 
the presence of an additional (but smaller) population of embed- 
ded YSOs which are not detected in the CCCP X-ray data. 



3. Color-Magnitude Diagrams for selected clusters 
in the CNC 

The stars in the CNC are distributed in a number of clusters 
(most notably the prominent optical clusters Tr 14, 15, and 16) 
and a widely dispersed population of young stars. In this sec- 
tion we consider the properties of the stellar populations in these 
different parts of the CNC. In order to avoid the background 
contamination problem, we restrict our analysis to the X-ray se- 
lected members, using the results of the clustering analysis of 
IFeigelson etal] d201 if . 



In our analysis of the CMDs we use the isochrones de- 
rived from the pre-main sequence models of lBaraffe et al.l ( 1998b 
for th e mass range 0.02 to 0.5 M Q and those of Siess et all 
(l2000h for the mass range 0.5 to 7 M . For the more massive 
stars, which are on the main sequence or already in their post- 
main sequence phases, we used the Geneva stellar models from 
iLeieune & Schaererl (12001 1) (model iso-c020-0650) for the mass 
range 7 to 70 M Q . 

Before we consider the individual CMDs, we want to point 
out that there are fundamental limitations to the accuracy of 
stellar mass- and age-determinations from CMDs (or HRDs); 
different studies of the same stellar population can sometimes 
produce sev erely discrepant results (see discussion in iHartmannl 
1200 ll 120031) . especially with respect to the presence or absence 
of age spreads. Factors such as differential extinction, photo- 
metric variability, overestimates of the stellar luminosity due 
to unresolved binary companions, and the effects of accretiorQ 



can ca use substantial scatter in the diagram (see, e.g., d iscus- 
sion in IPreibisch & Zinneckedfl999t ISlesnick et al.ll2008l) . The 
observed apparent luminosity spreads in the d ata are easily mis- 
interp reted as age spreads (see discussion in iHillenbrand et al.l 
120081) . It is therefore important to keep in mind that any observed 
luminosity spread in a CMD (or HRD) is always only an upper 
limit to a possible age spread. Underestimating the uncertain- 
ties in the observational data can easily lead to biased results^. 
We will therefore not attempt to determine ages and masses for 
individual stars, but restrict our interpretation of the CMDs to 
estimates of the typical ages for total population of young stars 
in each cluster; in a large enough sample, most of the above- 
mentioned uncertainties are expected to cancel out statistically. 
With these caveats in mind we now take a look at the CMDs for 
individual parts of the CNC. 



3.1. The cluster Trumpler 16 

The loose cluster Tr 16 is located in the center of the Carina 
Nebula and includes the optically dominant star r\ Car as well 
as the majority of the O-type stars in the CNC. All age esti- 
mates for Tr 16 reported so far in th e literature concern o nly the 
high- and intermediate mass stars. iMassev et al.l (1200 ll) found 
that the high-mass population of Tr 16 seems to be co-eval 
and report an age of 1.4 Myr (note that they assumed a dis- 
tance of 3.1 kpc rather than the value of 2.3 kpc we use here). 
iDeGioia-Eastwood et al. I J2001) claimed that the intermediate- 
mass stars in Tr 16 have been forming continuously over the 
last 10 My r f whereas t he hig h-mass stars formed within the 
Dias et ail (2002) list an age of 5 Myr, whereas 



last 3 Myr . 

iTapia et all (120031) suggested continuous star formation since 



6 Myr until now. These claims of large age spreads may be 
overestimates caused by the considerable differential reddening 
across the areas of these loose clusters. The most reliable age 
estimate seems to be the value o f ~ 3 - 4 Myr that was sug- 
gested bv lSmith & Brooksl d2008l) based on detailed stellar evo- 
lution models for r\ Car, the presence of 03 main sequence stars, 
and the modeling of the main-sequence turn-off in optical color- 
magnitude diagrams. 

The CMD of Tr 16 shown in Fig. g] is based on the 449 
stars with complete /- and H- and K s -band photometry among 
the 529 X-ray sources that are classified as Carina members 
and associated to one of the sub-clusters in the Tr 16 area by 
IFeigelson et alJd201 lb . Concerning the age of the low-mass stel- 
lar population, we note that our CMD of the X-ray selected 
sources shows a conspicuous group of stars with J ~ 12-14 
apparently aligned with the 3-4 Myr isochrones. The bulk of the 
low-mass stars is also seen at locations consistent with ages of 
« 3 - 4 Myr, and thus we consider this to be the most likely age 
for the low-mass stellar population in Tr 16. This value is consis- 
tent with the most likely age of the high-mass members, suggest- 
ing that high- and low-mass stars have formed at the same time. 
We note that a more comprehensiv e study of the X-ra y selected 
population of Tr 16 can be found in lWolket al.l (1201 lh . 



2 iBaraffe et alj d2009l) showed that stars with intermittently variable 
accretion rates can produce a large spread of luminosities such that a 
population of stars with identical ages of a few Myr may be wrongly 
interpreted to have an age spread of as much as ~ 10 Myr. 



3 To mention just t wo examples we ref er to the study of the young 
cluster NGC 3293 by Ba ume et all d2003l) . in which a careful analysis 
led to a corre ction of previous cla ims about the stellar ages, and to the 
new results of ICurrie et al.l (1201 Oh for h and x Persei, that clearly refute 
earlier claims of different ages and large age spreads for the double 
cluster. 
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3.2. The cluster Trumpler 14 

Tr 14 is the second most massive cluster in the Carina Nebula 
and contains 10 known O-type stars. Its spatial configura- 
tion is considerably more compact than Tr 16. The published 
age estimates for Tr 14 show re markable discrepancies: While 
iDeGioia-Eastwood et alJ d200l|) claim c ontinu ous star forma- 
tion ov er the last 10 My r, lDiasetal.1 (120021) list an age of 
2 Myr. iTapia et all d2003l) suggested a scenario of continuous 
star formation lastin g since 5 Myr until less than 1 Myr ago. 
lAscenso et al 1 (120071) claimed a more or less continuous age dis- 
tribution between zero and ~ 5 Myr, and suggested three peaks 
in the age distr ibution at values of 0-0.3 Myr, 1.4-2.5 Myr, 
and 3.2-5 Mvr. ISana et al.1 d20 1 Oh analyzed new adaptive-optics 
images of the cluster core (radius = 2') and derived a very young 
cluster age of only ~ 0.3-0.5 Myr. 

Our CMD of Tr 14 (Fig. |4]i is based on the membership to 
the rat her extended X-ray cluster as revealed by IFeigelson et al.1 
(1201 11) . It contains 1378 Chandra sources, for 1219 of which we 
have complete NIR photometry. The CMD suggests a somewhat 
broader age distribution than for Tr 16. Some objects lie between 
the 3 Myr and 10 Myr isochrones, but we note that most of these 
apparently older objects are located in the outer parts of the clus- 
ter. If we restrict the sample to the R — 2' cluster core, most 
of the apparently older objects are removed. The tendency for 
younger stel lar ages in the clu ster center agrees with the recent 
results from IS ana et al. The bulk of the stars in the clus- 

ter center seem to be quite young, < 3 Myr. Since this value is 
in reasonable agr eement with the ag e estimate for the high-mass 
stars in Tr 14 by iDias et al. 1 (120021) . we conclude that low- and 
high-mass stars in Tr 14 seem to be coeval. Many stars in the 
outer parts of the cluster seem to be somewhat older; however, it 
remains unclear whether these stars in the outer regions are ac- 
tually members of the cluster Tr 14, or whether they may be part 
of the widely distributed population of (partly older) stars in the 
CNC. 



3.3. The cluster Trumpler 15 

The cluster Tr 15 is located in the northern part of the CNC and 
contains 6 O-type stars. It is though t to be several Myr older than 
Tr 16 and Tr 14. [Dias et al.1 (120021) list an age of 8 Myr, while 
ITapia et al.1 d2003f claim a range of ages between ~ 4 Myr and 
~ 30 Myr, with a median value of about 8 Myr for the high-mass 
stars in Tr 15. 

. 4]i is again based on the membership analysis 
(201 ll). It contains 481 Chandra sources, for 



Our CMD (Fiji 
of IFeigelson et al.l 
436 of which we have NIR photometry. The upper part of the 
CMD is dominated by a group of stars lying between the 3 Myr 
and the 10 Myr isochrones; together with the clear lack of bright 
objects to the right of the 3 Myr isochrone, this suggests an age 
of ~ 5 - 8 Myr for the low-mass stars. This agrees with the age 
derived for the high-mass stars in Tr 15, suggesting that high- 
and low-mass stars are co-eval, and clearly showing that Tr 1 5 is 
older than Tr 16 and Tr 14. Finally, we note that a comprehensive 
study of the X-ray selected population of Tr 15 can be found in 
IWang et al.ld201ll) . 



3.4. The "Treasure Chest" cluster 

The most prominent of the embedded clusters in the Southern 
Pillars region is the "Treasure Chest", which is thought to be 
very young (< 1 Myr; see Smith et al. 2005). The most mas- 
sive member is a 09.5V star. The X-ray clustering analysis of 



IFeigelson et al ] (1201 lh associated 96 Chandra sources to the 
Treasure Chest. The 78 sources for which we have NIR photom- 
etry are plotted in Fig. |4] The typical extinction of the member 
stars, Ay > 5 mag, is considerably larger than in the other re- 
gions. The CMD of the X-ray selected stars is consistent with a 
very young age, since nearly all stars lie well above the 1 Myr 
isochrone. 



3.5. The cluster Collinder 232 

The clustering Cr 232 is located about 3' to the east of Tr 14 and 
defined by the early type stars HD 93250 (03.5V), HD 303311 
(05), and HD 93268 (A 2). Until rec e ntly, t he reality of a stel- 
lar cluster was not clear. ITapia et akl (|2003) noted that Cr 232 
may be just a random coincidence of these three relatively bright 
stars. However, the spatial distribution of the X-ray detected 
young stars in the CCCP clearly shows a sign ificant density 
peak in the Cr 232 region (IFeigelson et al.ll201 ll) . There are 70 
Chandra sources associated with this peak, for 67 of which we 
have J+H+K s photometry. The CMD is consistent with the as- 
sumption that the majority of the stellar population in this region 
is very young, < 3 Myr. 

Our HAWK-I images reveal several very red, deeply embed- 
ded infrared sources in Cr 232, including the pro minent circum- 
stellar disk object discussed in more detail in iPreibisch et al.l 
(2011c). This detection of embedded infrared sources, and the 
jets discussed in Sect. [7] clearly shows that Cr 232 actually is 
a cluster of very young stars and with ongoing star formation 
activity. We note that the IRAS and MSX data also led to the de- 
tectio n of deeply embedded luminous infrared sources in Cr 232 
(e.g., iMottram et al.l 120071) . some of which may be relatively 
massive YSOs. 



3.6. The cluster Bochum 1 1 

Bochum 1 1 (see Fig. |2| is a loose open cluster in the southeast- 
ern part of the CNC (i.e. in the South Pillars region) and contains 
5 O-type stars. From optical color-magnitude dia grams, age es- 
timates of < 3 M yr dFitzgerald & Mehtdfl987h and ~ 6 Myr 
dDias et al 1 120021) for the high-mass stars have been published. 

The X-ray clustering analysis associated 136 Chandra 
sources to Bochum 11. The CMD of the 117 X-ray selected stars 
with available NIR photometry (Fig.|4]i shows a group of stars ly- 
ing between the 3 Myr and the 10 Myr isochrones. This suggests 
that Bo 11 is somewhat older than Tr 14 and Tr 16, probably 
around 5 Myr. This agrees re asonably well wit h the age estimate 
for the high-mass stars from lDias et al.l(l2002l) . 



3.7. The "widely distributed" population 



According to the clustering analysis of IFeigelson et al 1 (120111) . 
5185 of the 10714 X-ray detected Carina members are not as- 
sociated with any cluster and thus constitute a widely dispersed 
population of young stars. Note that due to the spatial sensitivity 
variations (caused by the off-axis mirror vignetting and degrada- 
tion of the point spread function), the X-ray detection complete- 
ness for this dispersed population is somewhat lower than for the 
previously known clusters (that were preferentially observed at 
small on-axis angles in the CCCP survey mosaic). This implies 
that the relative size of this dispersed population (compared to 
the X-ray detected populations of the clusters) is somewhat un- 
derestimated. 
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Widely Distributed Population 




3.0 



Fig. 5. Color-magnitude diagram for the X-ray selected Carina 
members in the "widely distributed population" within the 
HAWK-I survey area. The meaning of the symbols is as in Fig. [4] 



For 1412 of these objects in the HAWK-I survey area we 
have complete /- and H- and X'j-band photometry. The CMD 
of these distributed stars (Fig. [5} is clearly different from the 
CMDs of the clusters Tr 14, 15, and 16, because it shows a much 
broader spread of colors. The considerable number of objects to 
the left of the 3 Myr isochrone and several objects near or to the 
left of the 10 Myr isochrone suggest a rather broad age distri- 
bution. We also note that the relative number of bright objects 
is considerably smaller than seen in the CMDs of the clusters; 
this suggests that the fraction of high-mass stars is smaller in the 
distributed population compared to the clusters. 

4. NIR excess fractions and ages of the different 
clusters in the CNC 

The fraction of X-ray selected stars w i th JHK excess emission 
has been determined in lPreibisch et al.l(l201 lbl) to be (9.7±0.8) % 
in Tr 14, (2.1+0.7)% in Tr 15, (6.9+1.2)% in Tr 16, and 
(32.1 ±5.3)% in the Treasure Chest. For Bochum 11 we find 
that 10 of the 117 sources with NIR photometry have NIR ex- 
cesses, giving an excess fraction of (8.5+2.6) %. Here we con- 
sider the relation between these excess fractions and the ages of 
the individual clusters derived from the analysis of the CMDs 
as described above. Figure [6] shows our results for the clusters 
in the CNC and compares it to similarly determined JHK excess 
fractions in other galactic clusters. 

It is obvious that the NIR excess fractions for all the clusters 
in Carina are systematically and considerably lower (by a fac- 
tor of about two) than those for the other young galactic clusters 
of similar ages. We interprete this difference as a consequence 
of the very harsh conditions due to massive star feedback in 
the CNC. The other clusters used for comparison contain much 
smaller numbers of massive stars and are thus characterized by 
much more quiescent conditions. The large population of very 
massive stars in the CNC leads to a level of hydrogen ioniz- 
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Fig. 6. JHK excess fraction versus age for the clusters in 
the CNC. For comparison, we also sho w JHK excess frac - 
tions derived for NGC 2024 and Taurus (iHaisch et al.l I20001 
the Orion Nebula Cluster ( ONC: IWolk et al.l 120051) . IC 348 
( IHaisch et alj|2001l), et a Cha dLvo et al.ll2003l) . Ori OB la and 
lb dBriceflo et al.ll2005l) . For the case of NGC 661 1, we plot the 
spatial variation of excess fract ions from the center of the cluster 
to the outer parts as derived bv lGuarcello et al.l (12007b . 



ing flux that is about 15 times higher than in the Orion Nebula 
(ISmifh & Broo ks 2008); the resulting increased heating and pho- 
toevaporation of the circumstellar disks around low-mass stars in 
the C NC is expect ed to lead to a quicker dispersal of these disks 
(e.g.. lClarkdl2007l) . 

This interpretation is supported by results about the spa- 
tial variation of the NIR excess fraction in the young cluster 
NGC 66 1 1 , that also c ontains a considerable p opulation of mas- 
sive stars in the center. iGuarcello et"aT1 (120071) found that the ex- 
cess fraction of the X-ray detected low-mass stars is correlated 
with the distance from the massive stars and drops from 24% for 
the outer parts of the cluster to 8% in the center (i.e. close to 
the massive stars). This variation agrees well with the difference 
in the excess fraction between the CNC clusters and the other 
galactic clusters of similar ages. 



5. The size of the stellar population in the CNC 

Our infrared data of the X-ray selected CNC members provide 
information about the masses of these stars that can be used to 
infer properties of the stellar mass function in the CNC. For this, 
we first have to consider the X-ray detection com pleteness. The 
study of the CCCP data by iBroos et all d2011al) shows that the 
completeness limit for lightly obscured stars observed at small- 
to-moderate off-axis angles on the detector is Lx ~ 10 30 erg/s in 
the 0.5 - 8 keV band. Assuming that the young stars in the CNC 
follow the relation between stellar mass and X-ray luminosity es- 
tabl ished by the data from the Chandra Orion Ultradeep Project 
(see lPreibisch et ai1l2005l) . we can expect a detection complete- 
ness of > 80% for stars with M* > 1 M Q . This completeness 
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drops towards lower masses and should become 50% for 
M« < 0.5 M e . 

We can thus assume the sample of X-ray selected Carina 
members to be nearly complete for stellar masses M, > 1 M . In 
order to determine the observed number of stars above this mass 
limit, we used the CMD for the X-ray selected Carina members 
in the HAWK-I survey area. We find 3182 objects at CMD posi- 
tions corresponding to M, > 1 M© for an assumed age of 3 Myr 
(i.e. above the reddening vector originating from the predicted 
position of 3 Myr old 1 M stars). 

It is interesting to compare this number to an IMF extrapo- 
lation based on the known popula t ion of high-mass stars in the 
CNC. The compilation of lSmifhl d2006h lists 127 individually 
identified members with known spectral types, and 52 of these 
should have stellar masses M > 20 M Q according t o the calibra- 
tions of iMartins et al.1 (120051) . The recent study of iPovich et alJ 
(2011b) lead to the detection of further candidates for young 
high-mass stars in the CNC and suggests that the true num- 
ber of high-mass stars should be about ~ 50% larger than the 
known population. Assuming thus a number of « 78 stars with 
M > 20 M , we can use the model representation of the field 
star IMF bv lKroupal d2002h and find that the predicted total pop- 
ulation of stars with At, > 1 M is ~ 4400. Since 80% of the 
known high-mass stars are inside the HAWK-I survey area, we 
thus have to expect that ss 3500 stars with M, > 1 M G should be 
found within the HAWK-I survey area, if the IMF in the CNC 
follows the Kroupa field star IMF. 

This prediction agrees very well to the number of 3 1 82 X-ray 
detected stars with CMD locations corresponding to M* > 1 M . 
This agreement directly shows that the number of X-ray detected 
low-mass stars in the CNC is as high as expected from field IMF 
extrapolations of the high-mass stellar population. It is a direct 
confirmation that the IMF in the CNC is consistent with the field 
star IMF down to (at least) M, * 1 M . 

This conclusion is further supported by the finding that the 
shape of the K-band luminosity function of the X-ray selected 
Carina member s agrees well with th at derived for the Orion 
Nebula Cluster dPreibisch et al. 201 lb). This implies that, down 
to the X-ray detection limit around 0.5-1 M , the shape of the 
IMF in Carina is consistent with that in Orion (and thus the field 
IMF). 

These results directly show that there is clearly no deficit of 
low-mass stars in the CNC down to ~ 1 M & . This is important 
because the issue whether the IMF is universal or whether there 
are systematic IMF variations in different environments is still 
one of the most fundamental open questions of star-formation 
theory (see lBastian et al.l2010i and references therein). It was of- 
ten claimed that some (very) massive star forming regions have a 
truncated IMF, i.e. contain much smaller numbers of low-mass 
stars than expected from the field IMF. However, most of the 
more rece nt and sensitive studies of massive star forming regions 
(see, e.g.. lLiu et al.ll2009t lEspinoza et alJl2009l) found the num- 
bers of low-mass stars in agreement with the expectation from 
the "normal" field star IMF. Our result for the CNC confirms 
this and supports the assumption of a universal IMF (at least in 
our Galaxy). In consequence, this result also supports the no- 
tion that OB associations and very massive star clusters are the 
dominant formation sites for the galactic fie ld star population, as 
already suggested bv lMiller & Scalol (11978b . 

Based on this result, we can now proceed and make an es- 
timate of the total stellar population (i.e. down to * 0.1 M Q ) of 
the CNC by assuming that the X-ray luminosity function in the 
CNC is similar to that in Orion (as sugg ested by the data). The 
study of the CCCP source statistics by iFeigelson et al.l d201 ll) 



found that the 3220 bright X-ray sources with a photon flux of 
\og(F) > -5.9 photons/sec/cm 2 in the 0.5-8 keV band consti- 
tute a spatially complete sample. If we transform this flux limit 
according to the difference in distances (i.e. 415 pc for the ONC 
versus 2.3 kpc for the CNC), the comparison to the distribution 
of photon fluxes in th e Chandra Orion Ultradeep Project sample 
dGetman et alJl2005l) yields an estimate of ~ 43 727 stars in total 
for the CNC. This number agrees very well with the extrapola- 
tion of the field IMF based on the number of massive stars (see 
above): for 78 stars with M t > 20 M G , the Kroupa IMF predicts 
that there should be « 40 000 stars with M« > 0. 1 M B . 

Finally, we can multiply the extrapolated number of ~ 
43 727 stars with the mean stellar mass of 0.64 M (as valid for 
the Kroupa IMF in the 0.1 M Q to 100 M range), and find a total 
stellar mass of the inferred CNC population of about 28 000 M . 

6. Search for faint star clusters and results on the 
spatial distribution of the young stars 

The spatial distribution of the stars in the CNC contains impor- 
tant information about the structure, dynamics, and evolution of 
the region. A very important question in this context is whether 
the stars form in a clustered mode or in a dispersed mode. Until 
recently, only the high-mass part of the stellar population was 
well known. Most of the massive stars are in one of the open 
clusters (and this was the reason why the CNC was often con- 
sidered to be a "cluster of clusters"). However, according to 
the field IMF, most of the total stellar mass is in the low-mass 
stellar population^. Recent studies have shown that a significant 
fraction of all stars in the solar neighb orhood seems to form i n 
a non-clustered, disp ersed mode (e.g.. iGutermuth et al.ll2009h . 
iBressert et al.l (1201 Ol) suggested that the stellar surface densities 
of young stars in nearby star forming regions follow one sin- 
gle smooth distribution rather than two different discrete modes 
(i.e. clustered versus distributed). 

The spatial distribution and clustering properties of the 
X-ray detected YSOs in the CNC was studied in detail by 
IFeigelson et al.l (1201 ll) . They identified 20 principal clusters of 
X-ray stars (most of which correspond to known optical clusters 
in the CNC) and 3 1 small groups of X-ray stars outside the ma- 
jor clusters. Altogether, these clusters contain about half of the 
X-ray detected YSOs in the CNC. The other half of the X-ray 
detected YSO population seems to constitute a widely dispersed, 
but highly populous, distribution of more than 5000 X-ray stars. 
A similar result was found in the analysis of the Spitzer data: 
most of the ~ 900 identified YSO candidates are spread through- 
out the South Pillar region (ISmith et al.ll2010bl) . Although this 
Spitzer study led to the detection of eleven previously unknown 
clusters, the populations of these clusters are so small (at most 
< 35 YSOs per cluster) that only a small fraction of all YSOs are 
in one of these new clusters. While the analysis of the Chandra 
and Spitzer data yielded consistent results, the limited sensitiv- 
ity of both data sets may leave some faint clusters undetected. 
The X-ray data are known to be seriously incomplete for masses 
S5 0.5 M . The same is true for the Spitzer data, that have a com- 
pleteness limit of [3.6] ~ 13, corresponding approximately to 
M» = 1 M . This implies that clusterings consisting of only a 
few dozen low-mass stars could have easily been missed. Such 
clusterings should, however, be easily visible in the HAWK-I im- 
ages, if their spatial configuration is compact enough. We have 

4 According to the Kroupa IMF, 89% of the total stellar mass in the 
0. 1 - 100 M range is contained in stars with M, < 20 M Q , and 73% in 
stars with M, < 5 M G . 
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Table 2. Parameters of the newly detected clusters. We list the 
apparent radius R, the number of stars N,, the number of NIR 
excess objects N %eX c> the number of X-ray detected stars N,,x in 
each cluster, and the /^-band magnitude of the brightest star. 



Name 


R 


N, jV,,^,, 


N,, x 


Brightest star 


[72000] 


["] 






K s [mag] 


C 104459-593118 


11 


* 33 8 


2 


10.90 


C 104513 -595753 


6 


* 10 1 


4 


9.30 


C 104559-600515 


8 


-16 3 


2 


10.60 


C 104613 -595832 


7 


w 14 8 


2 


12.86 



therefore performed a detailed visual inspection of the HAWK-I 
images to search for yet unrecognized clusters. 

The very inhomogeneous cloud extinction causes strong spa- 
tial variations in the surface density of observed background ob- 
jects in the HAWK-I images; this limits the ability to recognize 
loose clusterings consisting of only a small number (S; 10) of 
stars. We therefore restricted the search to clusters that con- 
sist (apparently) of at least 10 stars and have a configuration 
dense enough for them to clearly stand out from their surround- 
ings. This search revealed only four likely clusters that were not 
known before. Images of these clusters are shown in Fig. Qand 
their properties are listed in Table |2] We note that the physical 
nature of these apparent clusters is not entirely clear; some of 
them might be just random superpositions of unrelated objects at 
different distances. The cluster C 104559-600515 is the most re- 
liable, because this group of stars appears to be embedded in the 
head of a prominent pillar in the South Pillars region. However, 
even in this case it remains unknown whether the apparent clus- 
ter is a gravitationally bound physical entity. 

Despite the uncertainties about the physical reality of these 
clusters, our search provides a clear result: it directly confirms 
that there is no significant number of previously undetected clus- 
ters. The HAWK-I images clearly show that most of those young 
stars in the complex that are not associated to one of the already 
known clusters are in a non-clustered, dispersed spatial configu- 
ration. This result strongly su pports the conclusions d rawn from 
the X-ray clustering study of iFeigelson et al 1 (1201 lh that about 
half of the total young stellar population is in a widely distributed 
spatial configuration. It also supports the picture of small-scale 
triggered star formation by radiative (and wind) feedback from 
the massive stars (see discussion in Sect. [8}. 

7. Search for protostellar jets 

Jets and outflows are important signposts for embedded proto- 
stars that often remain undetected at NIR wavelengths, and can 
reveal the latest generation of the currently forming stars . When 
the outflowing material interacts with the environment, shocks 
are generated where the flows impact with surrounding dense 
clouds. In these shocks, hydrogen molecules can be collision- 
ally excited into higher ro-vibrational levels and the subsequent 
decay of these excited states causes strong emission lines in the 
NIR wavelength range. The 2.12 //m v — 1-0 S(l) ro-vibrational 
emission line of molecular hydrogen is a very convenient and 
widely used tracer of these shocks (e.g ., iMcCaughrean et all 
Il994tismifh et alJ2007tlDavis et al.ll2008h . This was the motiva- 
tion for obtaining HAWK-I images through a narrow-band filter 
centered on the 2.12 /^m Ht line. 

It is, however, important to keep in mind that the detection of 
ro-vibrational line emission is not a proof of the presence of jet- 
induced shocks, because molecular hydrogen can also be excited 



in higher r o-vibrational levels by the mechanism of UV fluores- 
cence fe.g. JBlack & Dalgarnoll976h . Since the Carina Nebula is 
characterized by very high levels of ionizing photon fluxes, UV 
fluorescence will be a very important excitation mechanism for 
H2 ro-vibrational emission. In our HAWK-I data, a distinction 
between collisional (shock) excitation and UV fluorescence can 
be made by comparing the strength of the emission seen in the 
H2 filter images to that in the images obtained through another 
narrow-band filter centered on the Bry line. Collisionally (jet-) 
excited molecular hydro gen from jets is no t expected to show 
Bry emission (see e.g., iNisini efal] 120021; ISmifh et alj 120071; 
iGarcia Lopez et all 1201 Oh . whereas UV excited (i.e. irradiated) 
matter should also show Bry emission. 

7.1. H 2 jets in the HAWK-I images 

We performed a detailed visual inspection of the HAWK-I 
narrow-band images to search for features that are bright in 
the H2 filter and much fainter (or invisible) in the Bry fil- 
ter. This search revealed only six clear cases of Molecular 
Hydrogen Emission-Line Objects (MHOs). These objects have 
been included in the "Catalogue of Molecular Hydrogen 
Emission-Line Objects in Outflows from Young Stars" 
(http://www.jach.hawaii.edu/UKIRT/MHCat/), where they are 
listed as MHO 1605 to 1610. 

The two most prominent H2 jets, MHO 1607 and 
MHO 1608, are found in the north-western part of the cluster 
Cr 232. A comparison of the H2 image and the Bry image is 
shown in Fig. [8] The fact that these features are bright in the 
2.121 yum H2 line but not visible in the Bry line confirms the na- 
ture of collision (i.e. jet-shock) induced emission. 

MHO 1607 (J104428. 2-593245) is composed of (at least) 
five discernable knots that extend over about 6", approximately 
in east-west direction. The morphology seems to suggest that the 
jet is moving towards the west. The HAWK-I images do not re- 
veal an obvious candidate for the source of the flow. The second 
object, MHO 1608 (J104429. 1-593242), is more compact and 
less bright. It seems possible that these two MHOs belong to a 
single flow system. In that case, the total linear extent is about 
1 1" and corresponds to a physical length of (at least) 0.12 pc. 

Fig. [8] shows that the surface of the clouds surrounding 
Cr 232 emits in the H2 line and also in the Bry line, showing 
that this is fluorescent emission from the ionization front at the 
surface of the cloud. In addition to the presence of the de eply 
embedded YSOs in this cloud (see iPreibisch et all 1201 id) , the 
detection of these jets clearly shows that Cr 232 is a prominent 
site of actively ongoing star formation in the CNC. 

Figure|9]shows the other four MHOs detected in the HAWK- 
I images. MHO 1605 (J10435 1.4-593911) is located south of 
Tr 14. The jet is embedded in the south-eastern tip of an hour- 
glass shaped globule. Two patches of H2 emission are seen above 
and below the rim of the globule, and a faint point-like source 
is seen inside the globule and just between them. MHO 1606 
(J104422.7-593354) is located south-west of Cr 232. It consists 
of a diffuse knot with a weak 2" extension to the north-east. 
MHO 1610 (J 104608 .4-594524) is located on the northern rim 
of a large pillar in the Tr 16 region. It consists of two diffuse 
blobs separated by 1.5". 

MHO 1609 (J1045 10. 1-600229) is located in the middle of 
a small elongated cloud in the western parts of the South Pillars 
region. This is the only case of a jet found in the central region of 
a globule. The projected length of the flow system is « 10", cor- 
responding to a physical length of (at least) 0. 1 1 pc. We note that 
this cloud is the source of the five optical jet candidates HHc-4 to 
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HHc-8 described bv lSmith et"aD (1201 Oal) . but the MHO seems to 
be unrelated to these optical HH objects. On the other hand, for 
none of the optical HH object candidates do we see any molecu- 
lar hydrogen emission. The HAWK-I images show a faint diffuse 
point-like source that may be related to HHc-6, but since this ob- 
ject is similarly bright in the H2 and Bry images, what we see is 
not H2 line emission but probably reflected light (perhaps the 
envelope of the jet-driving protostar). 



nario predicts that stars form very close to the irradiated edge of 
the clouds, where the radiative compression leads to collapse and 
triggers star formation. In the alternative scenario, if one would 
assume spontaneous star formation taking place at random loca- 
tions within these clouds, most protostars should be located in 
the inner parts of the clouds and we would thus expect to see 
more cases of MHOs within the clouds, which is not the case. 



7.2. Search for H 2 emission from the optical HH objects 



The HST Ha imaging survey of ISmith et al.l d2010al) revealed 
39 HH jets and jet candidates from YSOs that are embedded in 
dense globules. 13 of these HH jets and 14 jet candidates are lo- 
cated in the area covered by our HAWK-I narrow-band images. 
Our detailed inspection of the surroundings of the optical jets in 
the HAWK-I images showed that several of the optical HH jets 
can also be seen in the NIR narrow-band images, but in all cases 
the apparent brightness in the Bry line filter is similar to that in 
the H2 line filter. This implies that the emission we see in the H2 
filter image is not shock-induced but either related to UV fluo- 
rescence or just reflected continuum light. This is most likely a 
consequence of the fact that the optical HH jets seen by HST are 
irradiated atomic jets that have left the globules and expand in 
the surrounding diffuse atomic interstellar medium. This implies 
that no molecular hydrogen should be present in the immediate 
surroundings of the HH jets, and thus no H2 line emission is 
expected. 

However, one might expect to see H2 line emission from 
those parts of the jet flow that are still within the dense globules 
(i.e. the part of the path before the jets emerge into the surround- 
ing diffuse atomic medium). Therefore we performed a detailed 
inspection of those parts of the globules from which the optical 
HH jets emerge. If the jet-driving protostars were located deeply 
inside these globules, one would expect to see H2 line emission 
within the globule that should be aligned with the optical HH 
emission outside the surface of the globule. However, in no case 
do we find detectable H2 line emission at such locations. 

This null result is in agreement with the fact that none of 
the optical HH jets exhibited detectable excess 4.5 fim emission 
(which is also related to ro-vibrational line s from mole c ular hy - 
drogen) in the Spitzer imaged analyzed bv lSmith et al.1 (l2010bl) . 
This has an interesting implication: it shows that the driving 
sources of these jets are located very close to the edge of the 
globules, and not in the center of the globules (because then we 
would expect to see MHOs). For a protostar located at the edge 
(or at the tip) of a globule, the jets will immediately move into 
the diffuse atomic medium, and no MHOs are thus expected. 

7.3. Implications on the star formation mechanism 

We find only a very small number of six MHOs, i.e. signposts of 
protostellar jets, in the molecular clouds in the CNC. The clear 
lack of H2 jet detections at the positions of the numerous optical 
HH jets suggests that the driving sources are located very close 
to the surface of the globules, and the jets therefore move in 
the atomic gas, not through the molecular gas in the globules. 
This implies that almost all jet-driving protostars (i.e. currently 
forming stars) are located very close to the irradiated edges of 
the clouds, and not in the central regions of the clouds. 

This spatial configuration strongly supports the scenario that 
star formation is currently trigger ed by the advancing ionization 
front in the irradiated clouds (see lSmith et al. 2010b|. This sce- 



8. Summary and Conclusions 

8. 1 . Summary of the main results 

The analysis of the infrared properties of the young stars in the 
CNC provides important information about the stellar ages and 
masses, as well as the mass function and the total size of the 
stellar population. The ages estimated for the low-mass popula- 
tions in the young clusters within the CNC are consistent with 
previous age determinations for the massive cluster members, 
suggesting that the high- and low-mass stars have the same age, 
i.e. have formed together at the same time. 

The number of X-ray detected stars with CMD positions cor- 
responding to stellar masses of M* > 1 M Q is consistent with 
an extrapolation of the field star IMF based on the number of 
high-mass (M > 20 M Q ) stars. This suggests that the IMF in 
the CNC is consistent with the field star IMF (down to at least 
1 M G ). The extrapolation of the X-ray detected population down 
to 0.1 M Q suggests a total population of « 43 730 stars with an 
integrated mass of about 28 000 M . We note that this extrap- 
olation only considers the star numbers in the X-ray detected, 
lightly obscured population of young stars in th e CNC. The 
analysi s of deep Spitzer images of the CNC by iPovich et alJ 
(1201 lal) revealed a population of 1439 YSO candidates with 
strong mid-infrared excesses that is thought to be dominated by 
intermediate-mass (2 . . . 10 M ) stars. Most (72%) of these are 
not detected in the Chandra X-ray data, probably due to their 
strong obscuration. The extrapolations of IPovich et al.l (1201 lab 
suggests a total population of ~ 14 300 obscured objects, that 
have to be added to the X-ray selected population. This raises 
the total population of the CNC to » 58 000 stars, and the in- 
tegrated stellar mass to » 37 000 M Q . These numbers^ clearly 
show that the CNC constitutes one of the most massive galactic 
clusters/associations, probably topping the often quoted galac- 
tic "starburst temp lates" NGC 3603 (M» )to t ~ 15 000M o ; see 
iRochau et a"Hl2010h and th e Arches cluster (M«, tot ~ 20 000 M Q ; 
see lEspinoza e t al. 2009]), and being similar t o Westerlund 1 
(Af» >tot ~ 49 000M o ; see iGennaroet al.ll201ll) . The CNC can 
thus be regarded as the most nearby galactic cluster on the verge 
of extragalactic starburst clusters. 

The presence of several deeply embedded young stellar ob- 
jects and the detection of jets show that the cloud associated to 
the cluster Cr 232 is a center of ongoing star formation activity. 

The fact that we find only four new small clusters in the 
HAWK-I images confirm and extend previous results showing 
that about half of the total stellar population in the CNC is in a 
non-clustered, widely dispersed spatial mode. The presence of 
this large distributed stellar population suggests that the CNC 
should be considered as an unbound OB association rather than 
a cluster of clusters. Since most of the existing clusters are quite 



5 As a historical note, it is interesting to compare this numb er to older 
estima tes based on optical observations only: for example, iFeinsteir] 
i ll 9951) estimated the total stellar mass of the CNC to be 3596 M e . The 
progress in observations has lead to an increase by more than a factor 
of 10 within just about 15 years! 
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loose and will probably disperse within a few 10 Myrs, the frac- 
tion of widely distributed, non-clustered stars will even increase 
with time. 

The very small number of H2 jets we see in the clouds and 
the complete lack of molecular hydrogen emission related to the 
numerous known optical HH jets show that most of the very re- 
cently formed protostars must be located very close to the sur- 
face of the globules (not inside), suggesting a triggered star for- 
mation mechanism. 



8.2. Interesting aspects of the current star formation activity 
in the CNC 

The UV radiation and winds of the very hot and luminous early 
O-type and WR stars in the central area of the CNC affect the 
surrounding clouds strongly. While much of the cloud mass 
was (and still is) dissolved and streams away in the expanding 
super-bubble, the remaining molecular cloud material is now 
highly fragmented into numerous dense pillars. The total (dust 
+ gas) mass of the dense clouds in the CNC is still quite large 
(~ 60000M o ), but these cl ouds fill only a small fraction of 
the volume of the CNC (see iPreibisch et ai1l201 lab . The radia- 
tive compression of these clouds currently leads to triggered star 
formation. The spatial distribution of the YSOs detected in the 
Spitzer survey (ISmith et al.l2010bh suggests that the induced for- 
mation of stars happens predominantly near the most strongly ir- 
radiated tips of the clouds. This is consistent with recent results 
from numerical simulations of the evolution of irradiated clouds 
(e.g.. lGritschneder et al.ll2QToh . An implication of this triggered 
star formation mechanism is that the cloud material surrounding 
the newly formed protostars is simultaneously dissolved by the 
advancing ionization front. As their natal surrounding clouds and 
protostellar envelopes are dispersed very quickly, these newly 
formed stars are immediately exposed to the very harsh radia- 
tion field created by the luminous O-type stars, and therefore 
loose their circumstellar material very quickly. The transition 
from an embedded protostellar object to a revealed disk-less star 
thus happens faster than in star forming regions with lower levels 
of feedback. This can explain the very small fraction of young 
stars with NIR excesses as a tracer of circumstellar disks. 

Several aspects of this process of triggered formation of a 
new stellar generation in the CNC appear particularly interesting 
and are worth being discussed in more detail. First, the triggering 
process has produced many sparse groups of young stars instead 
of a few rich clusters; it is thus a qualitatively different mode 
of star formation than the process that formed the earlier gen- 
eration of stars in the much more populous clusters Tr 16 (and 
Tr 14). The relatively small sizes of the newly formed cluster- 
ings of young stars suggests that the currently ongoing radiative 
triggering of star formation is a very local, small-scale process. 
This is probably related to the highly fragmented structure of the 
now existing dense clo uds (as can be seen in the sub-mm map in 
IPreibisch et aT1l2011al) . The compression by the ionization front 
drives theses rather small individual clouds at various locations 
in the nebula into collapse, each of which then gives birth to a 
small stellar group. 

The second aspect is that the masses of the stars in the trig- 
gered new generation seem to be restricted to < 20 M Q , i.e. to 
much lower values than the very massive stars in the older, trig- 
gering generation (M ^ 1OOM for 77 Car and the 03 and WR 
stars in Tr 16 and Tr 14). The infrared and X-ray surveys re- 
vealed several candidates for intermediate-mass and moderately 
high-mass YSOs, but no trace of any very massive (M > 50 M G ) 
protostar has yet been found. This is again supported by the re- 



sults of the sub-mm st udy of the structure of the dense clouds 
(IPreibisch et al.ll201 lal) . where it was found that nearly all dense 
clouds in the CNC have masses < 1000 M Q ; according to the em- 
pirical relation between cloud mass and maximum stellar mass, 
these clouds are expected to yield maximum stellar masses of no 
more than ^15 M Q . 

The third aspect concerns the size of the triggered new stel- 
lar population. The X-ray survey showed that the widely dis- 
tributed population consists of a similar number (~ 5000 X-ray 
detected objects) of stars as the populations of the triggering 
populations in Tr 16 + Tr 14. Considering the spatial incom- 
pleteness effect discussed in Sec. 3.7, the dispersed population 
is probably even somewhat larger than the cluster populations. 
Evidence for a second generation of star formation triggered by 
radiative and wind-feedback from the massive stars in a first 
generation has been found in several other massive star form- 
ing regions, but in most cases the total size of the second stellar 
generation is considerably smaller than that of the first genera- 
tion (see. e.g.JStanke et al.l2002r.lReach et alJ2004tlLinskv et all 
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120071: IWang et al.ll2010h . The remarkably large size of the latest 
stellar generation in the CNC may be related to the particularly 
high level of massive star feedback in the CNC, that is perhaps 
more effective in triggering star formation than the much lower 
feedback levels in regions with smaller populations of very mas- 
sive stars in the older generation. 

The combination of the deep HAWK-I NIR images with the 
results of the X-ray survey has provided us with a very substan- 
tial amount of crucial new information about the stellar content, 
the star formation history, and the ongoing star formation pro- 
cess in the CNC. Future studies will combine these observa- 
tional data with detailed numerical simulations of how molecular 
clouds e volve under the influence of strong massive star feed- 
back (see lGritschneder et al.ll2010b . Such a comparison can pro- 
vide new and detailed insights into fundamental processes such 
as the disruption of molecular clouds by massive stars, the origin 
of the observed complex pillar-like structures at the interfaces 
between molecular clouds and H II regions, the effect of stellar 
feedback on molecular cloud dynamics and turbulence, and how 
ionizing radiation and stellar winds trigger the formation of a 
new generation of stars. 
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Fig. 4. Color-magnitude diagrams for the X-ray selected Carina members (crosses) in the different clusters in the CNC. Objects with 
infrared excesses are additionally marked by red open boxes. The solid green line shows the 3 Myr isochrone, the red and the blue 
lines show the low-mass (M < 7 M Q ) isochrones for ages of 1 Myr and 10 Myr, respectively. The large solid dots mark the positions 
of 1 M Q stars, the asterisks those of 0.075 M objects on these isochrones. The arrows show the extinction vector for Ay = 10 mag, 
and the dashed green line shows the 3 Myr isochrone redenned by that amount. The sequence of blue errorbars near the right edge 
indicates the typical range of magnitude dependent photometric uncertainties. In the CMD for Tr 14, the green crosses mark the 
X-ray sources in the central R = 2' cluster core. 
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Fig. 7. HAWK-I //-band images of the four newly detected clusters. 




Fig. 8. Images of the Cr 232 region obtained through the 2.121 fim H2 narrow -band filter (left) and the 2 . 1 66 //m Bry narrow -band 
filter (right). The location of the two molecular hydrogen jets found in this area, MHO 1607 (104428.2-593245) and MHO 1608 
(104429. 1-593242), are marked by the box and the circle in the upper righ t part of the images. Th e bright extended object just right 
of the center is the edge-on circumstellar disk object described in detail in lPreibisch et al.l (1201 id) . 
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Fig. 9. Images of the jets MHO 1005, 1606, 
1609, and 1610, obtained through the 2.121 yum 
H2 narrow-band filter (left) and the 2.166/jm 
Bry narrow-band filter (right). The location of 
the molecular hydrogen jets are marked by the 
arrows. 
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